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INTRODUCTION 


In  our  previous  report  (Rodriguez  et  al.,  1984),  we  presented  results 
from  our  diurnal  simulation  of  the  concentrations  of  the  sodium  species  based 
on  neutral  chemistry.  In  this  report,  we  extend  the  discussion  to  Include  ion 
chemistry  of  sodium  species. 

Positively  charged  ions  of  atomic  weight  23  have  been  observed  in  the  E 
region  during  several  rocket  flights  (see  Swlder,  1984,  and  references  there¬ 
in).  These  observations  have  been  interpreted  as  positively-charged  Na  ion 
(Na"*").  Sodium  ions  are  produced  by  photo-ionization  as  well  as  by  charge  ex¬ 
change  reactions  with  the  major  positive  ions  02^  and  NO'*’  (Swlder,  1969a). 
Dissociative  recombination  of  molecular  ions  produced  by  clustering  of  Na^ 
with  N2  and  CO2  (Richter  and  Sechrist,  1979)  ultimately  removes  Na^  Ions.  The 
clustering  rate  decreases  rapidly  with  altitude  resulting  in  Inefficient  re¬ 
moval  of  ions  at  high  altitudes.  Although  there  are  Insufficient  data  to  ob¬ 
tain  any  detailed  information  on  the  diurnal  and  seasonal  behavior  of  the 
ions,  they  can  nonetheless  affect  the  behaviour  of  the  neutral  sodium.  In 
particular,  the  departure  of  the  observed  scale  height  of  atomic  sodium  from 
the  scale  height  of  the  atmosphere  above  90  km  (Gibson  and  Sanford,  1970; 
Slmonich  et  al. ,  1979)  is  thought  to  be  controlled  by  the  amount  of  Na^  ions 
present  (Hanson  and  Donaldson,  1967). 

In  order  to  calculate  the  production  and  removal  rates  of  the  sodium 
ions,  one  needs  to  know  the  distributions  of  the  background  ions.  In  Sections 
2  and  3,  we  will  present  the  model  for  calculating  the  diurnal  behavior  of  the 
background  ions  in  the  D  and  lower  E  region.  Section  4  will  present  the  re¬ 
sults  of  the  calculated  diurnal  behavior  of  the  sodium  Ions.  The  effects  of 
the  ions  on  the  observations  of  the  scale  height  of  atomic  sodium  above  90  km 
and  other  Impact  on  sodium  will  be  discussed  in  Section  5. 
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2.  IONS  IN  THE  D  AND  LOWER  E  REGIONS 

Since  our  primary  interest  is  in  the  effect  of  metal  ions'  chemistry  on 
the  neutral  chemistry  around  90  km,  we  will  use  a  version  of  the  ion  chemistry 
appropriate  for  the  D  and  lower  E  regions.  In  the  lower  E  region,  production 
of  02^,  N2^  and  0'*'  by  x-ray  and  Lyman-6  dominates.  However,  the  N2^  and  0"'’ 
are  converted  rapidly  to  02^  by  charge  exchange.  Some  0^  is  converted  to  NO'*’, 
which  is  also  formed  by  the  charge  exchange  reaction  of  NO  with  02^.  In  the  D 
region,  production  of  NO^  by  Lyman-a  is  important,  making  this  ion  the  most 
important  ionic  constituent  near  80  km.  In  the  lower  D  region  (60-80  km), 
there  are  evidence  from  mass-spectrometrlc  measurements  that  hydrated  cluster 
ions  are  common  (Narclsl  and  Bailey,  1965;  Goldberg  and  Aikin,  1971;  Krankow- 
sky  et  al.,  1972).  Swlder  and  Narcisi  (1983)  suggested  that  carbon  dioxide 
may  be  an  important  agent  in  hastening  the  formation  of  these  clustered  ions. 

Negative  i  ms  are  produced  through  a  chain  of  reactions  initiated  by 
electron  attachment  to  O2 : 

e+02'*'M-*’02  +M.  (2—1) 

Although  observations  show  conflicting  results  making  the  identifications  of 
the  ions  difficult,  it  is  relatively  certain  that  the  major  ions  are  C03~  and 
NOa"  (Swlder  and  Narcisi,  1983).  For  the  problem  at  hand,  it  is  not  crucial 
that  the  identity  of  the  negative  ions  is  precisely  known. 

Positive  ions  are  destroyed  mainly  by  dissociative  recombination, 

AB^  +  e  -►  A  +  B, 
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where  AB'*’  denotes  a  molecular  Ion.  Rates  for  dissociative  recombination  are 
usually  very  fast  (~  10  ^  cm^  sec  1),  and  result  in  lifetimes  of  the  order  of 
minutes  for  positive  molecular  ions  near  90  km.  Removal  of  negative  ions  is 
effected  by  photodetachment  reactions 

B“  +  hv  B  +  e.  (2-3) 

As  the  concentrations  of  the  negative  ions  become  much  higher  below  80  km  at 
night,  positive-negative  ion  recombination 

A"*"  +  B~  +  products  (2-4) 

becomes  an  Important  loss  mechanism. 

A  simplified  model  of  the  reaction  schemes  for  positive  and  negative  ions 
appropriate  for  the  D  and  lower  E  region  are  given  in  Flgur  s  1  and  2  respec¬ 
tively.  The  photo-ionization  sources  Incorporated  in  our  model  include  the 
following; 

For  daytime  condition: 

a)t  X-rays ;  Absorption  and  ionization  cross  sections  of  0,  O2,  and  N2 
are  taken  from  Swider  (1969b),  for  wavelengths  between  1  and  lOOA. 
The  solar  fluxes  for  quiet  conditions  are  also  from  Swider  (1969b). 
b,c)  Lyman-6  ionization  of  O2  and  Lyman-g  ionization  of  NO:  Absorption 

cross  sections,  ionization  cross  sections,  and  fluxes  are  taken  from 
Swider  (1969a,  1983)  for  quiet  conditions, 
d)  Energetic  electrons:  Spectra  of  energetic  electrons  at  300  km  mea¬ 
sured  by  the  AE-C  satellite  have  been  presented  by  Torr  et  al. 

^Alphabet  refers  to  labels  on  Figure  3. 
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(1976).  These  authors  est  mated  an  ion  production  rate  of 
10  cm~3  at  for  nighttl  le  90  km  although  the  rate  for  daytime 
could  be  considerably  higher.  We  have  assumed  that  the  ion  produc¬ 
tion  profile  has  the  shape  of  the  theoretical  calculation  of  Berger 
et  al.  (1974;  Figure  6,  for  5  kev  e-folding  energy);  this  profile  is 
normalized  to  a  peak  production  rate  of  10  cm~3  at  90  km.  The 
above  total  production  rate  was  subdivided  into  90%  O2  and  10%  NO^ 
production  as  suggested  by  Swider  (1969a,  1983)  for  cosmic  rays. 
Because  of  the  large  uncertainties,  production  of  ions  from  ener¬ 
getic  electrons  is  not  included  in  most  of  our  calculations.  How¬ 
ever,  they  are  included  in  Figure  3  for  completeness. 

e)  Ionization  of  O2  (^A);  We  have  used  the  parameterization  suggested 
by  Paulsen  et  al.  (1972). 

f)  Cosmic  rays:  The  parameterization  given  by  Swider  (1983)  for  raid- 
latitude  condition  is  used. 


For  night  time  conditions: 

b,c)  Lyraan-8  and  Lyraan-g  scattered  by  the  geocorona:  The  scattered 

fluxes  as  a  function  of  altitude  and  solar  zenith  angle  are  taken 
from  Strobel  et  al.  (1974,  1980). 


d)  Energetic  electrons:  see  discussion  for  dayslde  ion  condition, 
f)  Cosmic  rays:  As  on  the  dayslde. 

A  constant  production  rate  of  10~3  cm~3  s~3  at  all  altitudes  due  to  galactic 
x-rays  (Torkar  and  Friedrich,  1983)  is  included.  The  production  of  ions  by 


starlight  (Strobel  et  al.,  1974,  1980)  has  been  ignored  since  it  is  only  im¬ 


portant  above  ~  110  km. 


A-4 


The  adopted  production  rates  for  ions  are  presented  for  noon  in  Figures 
3(a)  and  3(b)  and  3(c)  for  midnight.  These  calculations  assume  45°  latitude, 
equinox  conditions.  Ionization  by  x-rays  and  Lyman-g  are  dominant  above  an 
altitude  of  90  km  on  the  dayside.  The  nightside  sources  are  dominated  by  ion¬ 
ization  by  energetic  electrons  if  values  from  Torr  et  al.  (1976)  are  adopted. 
We  must  note,  however,  that  this  source  is  quite  variable  (Torr  et  al. ,  1976). 
The  profiles  for  ion  production  by  energetic  electrons  presented  in  Figure  3 
shjuld  thus  be  considered  as  upper  bounds.  Values  for  this  production  rates 
stniller  than  ionization  by  Lyraan-a  and  Lyman-6  cannot  be  ruled  out.  Uncer- 
t-i  ritles  also  exist  in  the  production  by  Lyman-a  and  Lyman-B  at  night,  since 
r*  •  'luxes  at  these  wavelengths  are  obtained  from  extrapolation  of  radiative 
;  r  uis^^-r  ra  1  fu lat ions  (Strobel  et  al. ,  1974,  1980). 

Tne  list  of  two-  and  three-body  reactions  included  in  the  model  is  given 
Ir  Table  I.  The  values  for  the  reaction  rates  are  taken  from  Swider  and  Nar- 
clsl  (  198  3),  Note  that  we  have  assumed  that  the  production  of  O2'*’  ions  is  due 
only  to  photo-ionization  processes.  The  additional  production  from  charge 
exchange  with  N2^  and  0^  is  Incorporated  into  the  same  term  (Swider,  1969a). 
The  production  of  NO'*’  is  from  both  photo-ionization  and  charge  exchange. 


3.  MODEL  FOR  THE  BACKGROUND  IONS 

The  basic  structure  of  our  ion  code  is  illustrated  by  the  flow  chart  in 
Figure  4.  This  code  was  adapted  from  our  mesopherlc  code  described  in  previ¬ 
ous  reports  (see,  e.g. ,  Rodriguez  et  al.,  1984).  The  local  concentrations  of 
the  electron,  the  positive  and  the  negative  ions,  are  governed  by  the  system 
of  equations: 


(3-1) 
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[e]  +  [NSUM]  =  [PSUM] 


(3-2) 


where  ,  and  [N^]  are  the  local  production  rate,  loss  frequency,  and  num¬ 

ber  density  of  the  1^^  ionic  species,  [e]  is  the  electron  density,  and  [NSUM] 
and  [PSUM]  are  the  sum  of  the  negative  and  positive  ions  concentrations,  re¬ 
spectively.  Equation  (3-1)  corresponds  to  the  local  approximation  of  the  con¬ 
tinuity  equation  for  each  ionic  species  whereas  (3-2)  serves  to  determine  [e] 
by  charge  conservation.  The  numerical  scheme  adopted  to  solve  the  system  can 
be  used  for  any  arbitrary  grouping  of  the  ionic  species  and  does  not  assume  a 
chain  structure  for  the  ion  chemistry. 

Equation  (3-1)  is  solved  at  a  series  of  discrete  time  steps  using  the 
fully  implicit  scheme  described  by  Sze  et  al.  (1980).  The  initial  values  for 
the  electron  and  positive  ion  densities  at  noon  are  taken  from  our  steady 
state  calculation.  At  other  time  steps,  the  Initial  values  are  taken  to  be 
the  electron  and  positive  ion  densities  calculated  at  the  previous  time  step. 
At  each  time  step,  the  equations  arc  solved  using  an  iteration  scheme  similar 
to  the  one  used  by  Swider  and  Foley  (1978).  The  values  for  [e]j+i  and 
[PSUM]j+^  for  the  j+1^^  iteration  aie  updated  via  (Swider  &  Foley,  1978) 

[e]j+^  =  {[e]j  +  [PSUM]j/(l  +  X)}/2, 

[PSUM]j^.j  =  {[e]j  +  [NSUM]j  +  [PSUM]j}/2  (3-3) 

with 

X  =  [NSUM]j/[e]j. 

The  process  is  repeated  until  equations  (3-3)  are  satisfied  to  some  specified 
accuracy.  The  solution  is  propagated  forward  using  time  steps  of  one  hour 
during  the  day,  two  hours  at  night,  and  30  minutes  around  sunrise  and  sunset. 
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After  24  hours,  the  solutions  are  tested  for  periodicity,  and  the  calculation 
is  repeated  if  periodicity  is  not  achieved  to  within  IZ.  Due  to  the  small 
time  constants  for  ionic  processes,  our  solutions  usually  become  periodic  in 
one  to  two  diurnal  cycles. 

The  neutral  background  atmosphere  is  obtained  from  the  diurnal  calcula¬ 
tion  described  in  previous  reports  (see,  e.g.,  Rodriguez  et  al.,  1984).  The 
effect  of  the  ion  chemistry  on  neutral  background  trace  gases  is  negligible 
under  most  conditions. 

Results  of  our  calculations  for  the  background  ions  and  electrons  at  45‘’N 
for  equinox  condition  are  presented  below.  The  diuinal  variation  of  the  elec¬ 
tron  density  is  presented  in  Figure  5  at  altitudes  of  86,  90,  94,  and  100  km. 
Similar  results  for  the  same  altitudes  are  shown  in  Figure  6  for  oj  and  NO^, 
which  are  the  main  positive  ion  constituents  in  this  region.  The  electron 
densities  exhibit  a  day-to-night  ratio  of  about  20  -  30  between  86  and  94  km. 
Increasing  to  about  50  at  100  km.  There  is  consideiable  variation  during  the 
daytime,  with  the  electron  density  increasing  by  factors  of  2  to  5  between 
sunrise  and  noon.  NO^  densities  are  larger  than  0^  densities  over  most  of 
this  region;  the  diurnal  behavior  of  NO^  is  thus  positively  correlated  with 
that  of  the  electrons.  Below  95  km,  0^  densities  exhibit  a  larger  day-to¬ 
night  ratio  (greater  than  200).  This  is  due  to  the  conversion  of  O2  to  NO'*’  at 
night  through  charge  exchange  with  NO. 

Our  results  during  the  day  are  in  good  agreement  with  diurnal  calcula¬ 
tions  presented  by  Torkar  and  Friedrich  (1983).  We  note  that  the  magnitude  of 
the  diurnal  variation,  however,  will  depend  on  the  model's  assumptions  for 
nightside  sources.  The  ionization  by  energetic  electrons  exhibits  great  vari¬ 
ability.  In  these  calculations,  we  neglect  ion  production  by  energetic  elec¬ 
trons.  The  Lyman-a  and  Lyman-6  sources  are  computed  from  radiative  transfer 
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models  which  become  uncertain  below  ~  100  km  (Strobel  et  al.,  1974,  1980). 

Our  results  should  be  Interpreted  only  as  the  general  diurnal  behavior  of 
electrons  and  ions.  More  detailed  information  on  the  nightside  ionization 
sources  at  specific  times  is  needed  to  pinpoint  the  diurnal  variation  to  a 
greater  accuracy. 

Profiles  for  electron  and  total  negative  ion  densities  for  day  and  night 
are  shown  in  Figure  7  for  the  same  conditions  as  in  Figures  5  and  6.  In  light 
of  the  possibly  large  variability  in  the  electron  flux,  we  also  include  the 
calculated  profiles  for  which  the  ionization  from  energetic  electrons  at  night 
is  included  (curves  e,  f).  The  figure  Illustrates  that  variability  of  this 
source  can  induce  changes  in  the  electron  density  of  more  than  factors  of  4 
above  90  km. 

The  distributions  of  negative  ions  for  daytime  and  nighttime  conditions 
are  given  in  Figure  8.  The  total  negative  ion  density  in  the  daytime  remains 
smaller  than  the  electron  density  down  to  about  60  km.  During  the  night,  the 
concentration  of  the  negative  ions  becomes  larger  in  the  absence  of  photode- 
tachraent  reaction  (2-4).  The  density  of  NO3  grows  to  more  than  10^  cm~3, 
reflecting  the  increase  in  concentrations  of  NO2  and  O3  which  favors  the 
transfer  to  NO3  via  O3  and  CO3.  Hydrated  negative  ions  are  not  included  in 
the  calculation  since  the  rates  for  the  reactions  involving  these  ions  are 
uncertain. 

Profiles  of  positive  ions  for  day-  and  nighttime  conditions  are  given  in 
Figures  9  and  10.  The  positive  ions  are  mostly  in  the  form  of  O2  and  NO"*” 
above  90  km.  Below  90  km,  various  hydrated  ions  are  formed  in  appreciable 
amounts.  Densities  for  these  ions  are  shown  in  Figure  10.  Formation  of  these 
hydrated  ions  is  Initiated  by  the  clustering  reactions  (Swider  and  Narcisl , 


(3-4) 


•  . 
B  a 

s 


HaO'*'  +  N2  -  m  HaO"^  •  N2  +  m, 

NO'*’  +  N2  +  N2  NO'*’  •  N2  +  N2, 
no"*"  +  CO2  +  N2  -*■  no"'"  •  CO2  +  N2, 
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followed  by  switching  reactions  of  the  form 

.  N2  +  CO2  -*•  a'*'  •  CO2  +  N2,  (3-5) 

A  •  CO2  +  H2O  >  A  •  H2O  +  C02« 

The  profiles  in  Figures  9  and  10  are  calculated  using  the  H2O  mixing  ratios  of 
Bevilarqua  et  al.  (1983),  which  are  up  to  an  order  of  magnitude  smaller  than 
those  of  Allen  et  al.  (1984)  above  80  km.  Variations  in  upper  mesospheric 
water  will  affect  the  calculated  densities  of  the  hydrated  positive  ions.  The 
results  in  Figures  9  and  10  reproduce  the  general  characteristics  of  calcula¬ 
tions  by  other  models  (Swider  and  Foley,  1973;  Swider,  1983;  Swider  and  Nar- 
cisl,  1983). 

4.  DIURNAL  MODEL  FOR  SODIUM  NEUTRALS  AND  IONS 

The  ionic  sodium  species  considered  are:  Na  ,  NaO  ,  Na02  ,  Na  *  N2, 

Na"*”  •  CO2 ,  and  Na^  •  H2O.  The  chemical  scheme  for  this  family  is  illustrated 
in  Figure  11.  This  scheme  combines  the  reactions  proposed  by  Murad  (1978)  and 
Richter  and  Sechrlst  (1979).  The  latter  paper  proposed  the  clustering  scheme 
that  recycles  the  Na^  back  to  neutral  Na. 

The  Na'*'  ion  is  produced  by 

Na  +  hv  Na^  +  e,  (4-1) 

Na  +  O2  ->■  Na"*"  +  O2  ,  (4-2) 

and 

Na  +  NO"^  +  Na'*'  +  NO.  (4-3) 
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•*-  •*. 


Molecular  ions  are  formed  via 


Na'*’  +  Nj  +  M  -»■  Na'*’*  Nj  +  M  (4-4) 

which  in  turn  yields  the  more  stable  Na  •  CO2  and  Na  •  H2O  through  the  switch¬ 
ing  reactions: 

Na"^*  Ng  +  CO2  +  Na‘‘‘«  COj  +  Nj  (4-5) 

and 

Na"*"*  CO2  +  HjO  -*■  Na'*’*  H2O  +  COj.  (4-6) 

The  cluster  ions  are  converted  back  to  neutral  atomic  sodium  by  dissociative 
recombination.  Reaction  (4-4)  is  the  rate  limiting  step  for  the  sequence  of 
reactions  (4-4)  through  (4-6). 

The  ion  reactions  shown  in  Figure  11  are  added  to  our  diurnal  model  for 
the  sodium  species  (Rodriguez  et  al.,  1984),  and  the  diurnal  behavior  and  al¬ 
titude  distributions  of  the  ionic  and  neutral  sodium  species  are  calculated  In 
a  self-consistent  manner.  In  these  calculations,  the  background  neutral  and 
ionic  constituents  are  not  recalculated,  since  they  are  not  affected  by  the 
sodium  chemistry.  Parameters  used  in  calculating  the  reaction  rates  in  our 
sodium  model  are  shown  in  Table  2,  with  the  rates  for  the  neutral  species  as 
documented  in  Rodriguez  et  al.  (1984)  and  those  for  the  ionic  species  taken 
from  Murad  (1978)  and  Richter  and  Sechrlst  (1979).  As  in  the  case  of  the 
neutral  chemistry,  most  of  the  ion  reactions  have  not  been  measured,  and  the 
values  proposed  are  deduced  from  measured  values  of  similar  reactions. 


The  diurnal  variation  of  Na^  is  shown  in  Figure  12  for  90,  94,  and  100 
km.  The  day-to-night  ratio  ranges  from  400  (90  km)  to  3  (100  km).  The  factor 
of  10  variation  for  the  altitude  range  90  to  110  km  suggested  by  the  data  as 
presented  by  Swider  (1984)  is  thus  consistent  with  our  results,  although  the 
large  changes  in  diurnal  variation  over  a  small  altl tude  range  preclude  draw¬ 
ing  any  definitive  conclusions  about  the  validity  of  our  chemical  scheme. 

Altitude  profiles  at  noon  for  Na^  and  Na  are  shown  in  Figure  13.  Atomic 
sodium  is  converted  into  Na"*"  at  high  altitudes  where  their  densities  are  com¬ 
parable.  As  we  shall  see  in  the  next  section,  processes  that  affect  the  con¬ 
centration  of  Na^  can  then  indirectly  alter  the  distributions  of  Na  at  these 
altitudes.  Densities  of  the  other  ionic  species  are  smaller  by  factors  of  100 
or  more.  In  particular,  the  [NaO"*"]/ [Na^l  value  is  ~  10”^  to  10~®,  in  agree¬ 
ment  with  the  estimations  of  Murad  (1978). 

5.  SODIUM  lOH  CHEMISTRY  AMD  THE  SCALE  HEIGHT  OF  THE 

NEUTRAL  ATOMIC  SODIUM  LAYER 

The  topside  scale  height  of  the  sodium  layer  has  been  observed  to  be  much 
smaller  than  that  of  the  background  atmosphere  (Gibson  and  Sandford,  1971; 
Megie  and  Blamont,  1977;  Slmonich  et  al.,  1979;  Granler  and  Megie,  1982),  de¬ 
creasing  to  about  2-4  km  at  100  km  altitude.  A  possible  explanation  for  this 
snail  scale  height  is  the  conversion  of  Na  to  Na^  by  photoionization  or  charge 
exchange  with  02^  (Hanson  and  Donaldson,  1967;  Gadsden,  1968;  Richter  and 
Stchrlst,  1979). 

This  sensitivity  of  the  scale  height  of  atomic  Na  ion  chemistry  is  illus 
trated  in  Figure  14  which  gives  the  calculated  profiles  of  Na  with  different 
assumptions  for  the  Na^  abundance.  As  discussed  in  Section  4,  reaction  (4-4) 
constitutes  the  rate-limiting  step  for  recombination  of  sodium  ions  and  is 
assumed  to  be  of  the  form  (Richter  and  Sechrist,  1979): 


A-11 


cm®  s~^. 


(5-1) 


Ks  =  K* 


A  value  of  K*  =  2  x  lO'^l  cm®  s”^  was  suggested  by  Richter  and  Sechrist  (1979) 
by  analogy  to  similar  reactions  involving  NO'*’.  Reaction  (4-4)  has  also  been 
measured  at  130°  K  in  a  helium  buffer  by  Perry  et  al.  (1980);  their  value  of 
2.6  X  10~3®  cm®  s~^  implies  a  value  of  K*  =  7  x  10”^^  cm®  s~^  if  the  tempera¬ 
ture  dependence  in  eq.  (5-1)  is  valid.  The  Na  profiles  in  Figure  14  are  cal¬ 
culated  using  different  values  for  K*.  The  noontime  profile  from  the  calcula¬ 
tion  without  ionic  species  (-o-)  is  compared  with  the  cases  calculated  with 

value  of  K*  from  Richter  and  Sechrist  (1979)  (K*  =  2  x  10~31  cm®  s''^,  _ )  and 

Perry  et  al.  (1980)  (K*  =  7  x  lO”^^  cm®  - ).  We  have  also  considered 

the  hypothetical  situation  with  K*  =  4  x  10~32  cm®  s“^  (IQiK).  As  expected, 
the  scale  height  of  neutral  sodium  increases  with  increasing  K*,  reflecting 
the  more  efficient  rate  for  reconversion  of  Na'*’  to  Na.  The  observed  small 
scale  heights  about  90  km  could  thus  be  explained  by  the  ion  chemistry  without 
having  to  invoke  transport  processes  (as  in  Hanson  and  Donaldson,  1967,  for 
example)  if  Kg  is  sufficiently  fast. 

Local  values  for  the  topside  scale  height  are  given  in  Table  3  for  the 
cases  considered  in  Figure  14.  A  comparison  with  observations  indicates  an 
upper  limit  of  2  x  10”31  cm®  s~^  for  K*  if  the  small  topside  scale  heights  aie 
to  be  explained  solely  by  ion-neutral  chemistry.  Such  small  values  of  K*  are 
not  likely,  however,  given  the  measurements  of  Perry  et  al.  (1980)  and  that 
the  faster  clustering  rate  are  expected  when  N2  is  the  buffer  gas  (Jegou  et 
al.,  1984;  Heimerl,  1976).  An  alternative  mechanism  to  reduce  the  topside 
scale  height  is  the  downward  transport  of  Na'*’  from  the  region  prior  to  its 
conversion  to  Na.  Vertical  ion  drifts  with  velocities  of  10  to  50  cm  s~^  ha'-e 
been  suggested  by  Hunten  (1981)  ani  Kirchoff  and  Batista  (1984)  and  could  be 
effective. 
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Ionic  species  could  acquire  a  vertical  velocity  through  electromagnetic 
interactions.  These  Include  electrostatic  forces  from  the  ionospheric  elec¬ 
tric  fields  (Kirchoff  and  Batista,  1984)  and  v  x  B  type  forces  from  the  mag¬ 
netic  field  because  of  the  motions  in  the  horizontal  plane  associated  with  the 
zonal  wind  (Jegou  et  al.,  1984)  and  tides  (Kirchoff  and  Batista,  1984).  Ion 
drift  could  provide  an  effective  way  for  the  vertical  transport  of  Na  species 
in  the  mesosphere,  particularly  when  a  significant  portion  of  the  sodium  spe¬ 
cies  is  in  ionic  form. 

A  detailed  calculation  of  these  effects  would  require  the  Incorporation 
of  the  above  parameters  in  our  model.  In  lieu  of  a  full  calculation,  we  exam¬ 
ined  the  response  of  the  topside  layer  to  typical  net  vertical  downward  velo¬ 
cities  for  ions  by  modifying  the  continuity  equation  for  the  NaX  species  to 
ha._  the  form 


d[([Na+]/(NaXj)  .  (NaX}  .  Wq] 
dz  ^  dz 


L, 


(5-2) 


where  is  the  flux  due  to  molecular  and  eddy  diffusion,  P  and  L  denote  the 
production  and  loss  mechanisms  for  NaX,  and  Wp  represents  the  vertical  ion 
drift  velocity.  The  relative  abundance  [Na^]/[NaX]  is  independent  of  the  as¬ 
sumed  drift  velocity  if  the  time  constant  for  the  vertical  drift  is  larger 
than  the  time  constant  for  the  removal  of  Na"*”  as  given  by  the  rate  limiting 
step  (4-4).  This  condition  is  valid  below  100  km,  where  the  time  constant  for 
the  vertical  drift  is  ~  10**  s,  while  the  time  constant  for  Na^  clustering  is 
~  103  s. 

Figure  15  shows  the  results  of  our  calculations  for  constant  drift  velo¬ 
cities  equal  to  0,  -5,  -10,  and  -20  cm/s  between  120  and  90  km.  The  small 
topside  scale  height  observed  (2-3  km,  Gibson  and  Sandford,  1981;  Clemesha  et 
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al. ,  1982)  can  be  modeled  using  dovmward  velocities  between  10  and  20  cm/s. 

The  effect  of  vertical  ion  drift  is  thus  an  important  contributing  factor  In 
controlling  the  topside  sodium  layer. 

6.  CONCLUDING  REMARKS 

We  have  presented  the  results  of  a  fully  diurnal  model  for  ion  species  in 
the  D  and  lower  E  regions.  When  coupled  with  our  previous  model  of  the  neu¬ 
tral  species,  our  present  model  allows  for  the  self-consistent  calculation  of 
background  neutral  and  ionic  species  between  60  and  120  km.  Results  of  our 
calculations  are  in  good  agreement  with  those  of  other  models. 

Production  of  Na"*"  occurs  primarily  through  photoionization  below  100  km, 
and  through  charge  exchange  with  02^  at  higher  altitudes.  The  main  loss 
mechanism  is  dissociative  recombination  of  sodium  cluster  ions.  A  self- 
consistent  calculation  thus  requires  knowledge  of  the  O2''’  ani  e  densities  in 
the  above  altitude  regimes.  Our  code  has  been  applied  to  the  calculation  of 
the  density  and  diurnal  variation  of  sodium  species  above  90  km.  The  calcu¬ 
lated  diurnal  variation  of  Na'*’  (factor  of  10-100  between  90  and  100  km)  and 
daytime  densities  (~  10^  cm '3)  are  in  good  agreement  with  rocket  data  (Swider, 
1984). 

We  have  also  lnvestiga:ed  the  effect  of  ion  chemistry  on  the  abundance  of 
neutral  sodium  above  90  km.  In  particular,  the  small  scale  height  (2-3  km) 
observed  in  the  topside  sodium  layer  seems  to  require  either  downward  ion 
drifts  velocities  of  order  10-20  cm  s~l  above  90  km  or  a  rate  of  clustering  of 
Na^  with  N2  a  factor  of  5  smaller  than  the  value  suggested  by  Richter  and 
Sechrist  (1979).  Further  elucidation  of  this  question  will  require  measure¬ 
ments  of  the  clustering  rate  at  mesospheric  temperatures  as  well  as  incorpora¬ 
tion  of  more  realistic  scenarios  for  ion  drifts  into  our  model. 
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Since  we  have  included  a  complete  treatment  of  the  ion  chemistry  in  the  D 


region,  the  model  presented  above  is  also  suited  to  studies  of  the  D  region 
ion  chemistry  and  its  possible  effect  on  neutral  constituents  during  disturbed 
conditions  (Swider  and  Keneshea,  1973;  Solomon  et  al . ,  1983).  The  model  also 
allows  for  each  extension  to  the  calculation  of  other  metallic  ion  species. 


Table  1 

List  of  Reactions  and  Reaction  Rates  for  the  Background  Ions 


TYpga 

PARI 

PAR2 

PAR3 

1 

7. 

l.OOE-28 

-3.00E+00 

O.OOE+00 

2 

1. 

l.OOE-09 

O.OOE+00 

O.OOE+00 

3 

1. 

4.50E-10 

O.OOE+00 

O.OOE+00 

4 

7. 

1.60E-30 

-3.20E+00 

O.OOE+00 

5 

1. 

3.00E-10 

O.OOE+00 

O.OOE+00 

6 

1. 

1.50E-09 

O.OOE+00 

O.OOE+00 

7 

1. 

l.OOE-09 

O.OOE+00 

O.OOE+00 

8 

1. 

2.00E-10 

O.OOE+00 

O.OOE+00 

9 

1. 

1.40E-09 

O.OOE+00 

O.OOE+00 

10 

7. 

1.40E-30 

-2.00E+00 

O.OOE+00 

1  1 

1. 

l.OuE-09 

O.OOE+00 

O.OOE+00 

12 

1. 

1 .OOE-09 

O.OOE+00 

O.OOE+00 

13 

7. 

5.00E-28 

-3.00E+00 

O.OOE+00 

14 

1 . 

l.OOE-09 

O.OOE+00 

O.OOE+00 

15 

7. 

3.50E-27 

-3.00E+00 

O.OOE+00 

16 

7. 

2.20E-27 

-3.00E+00 

O.OOE+00 

17 

7. 

2.30E-27 

-3.00E+00 

O.OOE+00 

18 

7. 

l.OOE-29 

-2.30E+00 

O.OOE+00 

19 

1. 

l.OOE-09 

U.OOE+00 

O.OOE+00 

20 

7. 

1 . lOE-27 

-4.70E+00 

O.OOE+00 

21 

7. 

1.60E-27 

-4.70E+00 

O.OOE+00 

22 

1. 

7.00E-11 

O.OOE+00 

O.OOE+OO 

23 

7. 

4.00E-07 

-1  .OOE+00 

O.OOE+00 

24 

7. 

2. lOE-07 

-6.30E-01 

O.OOE+OO 

25 

1. 

2.00E-06 

O.OOE+00 

O.OOE+OO 

26 

1. 

1 .50E-06 

O.OOE+00 

O.OOE+OO 

REACTIONS 

H'^.(H20)3  +  CO2  +  M  =  H‘^.(H20)3.C02  +  M,  7 
H'*’.(H20)3  .CO2  +  H2O  =  H'''.(H20)4  +  CO2,  1 
02'’'  +  NO  =  NO"*"  +  O2,  1 
02"*"  +  O2  +  M  =  04"^  +  M,  7 
04''‘  +  0  =  02''’  +  O3,  1 
04^  +  H2O  =  02^«H20  +  O2  j  1 
02'*‘.H20  +  H2O  =  H''’.H200H  +  O2 ,  1 

02''’.H20  +  H2O  =  H'^.H20  +  OH,  1 

H'’'.H200H  +  H2O  =  H'^.IH20)2  +  OH,  1 

H'^.H20  +  N2  +  N2  =  H'''.H20.N2  +  N2,  7 

H'^.H20.N2  +  CO2  =  H'^.H20.C02  +  N2,  1 
H'*’.H20.C02  +  H2O  =  H'^.(H20)2  +  CO2 ,  1 
H'^.(H20)2  +  CO2  +  M  =  H'^.(H20)2C02  +  M,  7 
H'^.(H20)2C02  +  H2O  =  H'*‘.(H20)3  +  CO2,  1 
H'^.H20  +  H2O  +  M  =  H'*'.(H20)2  +  M,  7 
H'''.(H20)2  +  H2O  +  M  =  H‘''.(H20)3  +  M,  7 

H'^.(H20)3  +  H2O  +  M  =  H''’.(H20)4  +  M,  7 

NO'^  +  CO2  +  N2  =  N0‘''.C02  +  N2,  7 
N0‘^.C02  +  H2O  =  CO2  +  N0'''.H20,  1 
N0'''.H20  +  H2O  +  N2  =  N0‘''.(H20)2  +  N2,  7 
N0'^.(H20)2  +  H2O  +  N2  =  N0'''.(H20)3  +  N2,  7 

N0’^.(H20)3  +  H2O  =  H'^.(H20)3  +  XX,  1 

NO"*"  +  E  =  N  +  0,  7 
02''’  +  E  =  0  +  0,  7 
04"*”  +  E  =  O2  +  O2  ,  1 
O2  ^  •  H2  0  +  E  =  O2  +  H2  0 ,  1 


The  parameterization  and  values  for  the  reaction  constants  are  taken  from 
Swider  and  Narclsi  (1984).  The  reaction  rates  are  calculated  as  follows: 


Type  1 
Type  7 
Type  8 


k  =  PARI 

k  =  PAR1*(300/T)"^^‘^^ 

k  =  PARl*(300/T)"PAR2*g^p(p^j^3/^^ 
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TYPE 

PARI 

PAR2 

PAR3 

27 

1. 

2.00E-06 

O.OOE+00 

O.OOE+OO 

28 

1. 

1.30E-06 

O.OOE+00 

O.OOE+OO 

29 

1. 

2.50E-06 

O.OOE+00 

O.OOE+OO 

30 

1. 

3.00E-06 

O.OOE+00 

O.OOE+OO 

31 

1. 

3.60E-06 

O.OOE+00 

O.OOE+OO 

32 

1. 

1.50E-06 

O.OOE+00 

O.OOE+OO 

33 

1. 

2.00E-06 

O.OOE+00 

O.OOE+OO 

34 

1. 

3.00E-06 

O.OOE+00 

O.OOE+OO 

35 

1. 

1.50E-06 

O.OOE+00 

O.OOE+OO 

36 

1. 

1.50E-06 

O.OOE+00 

O.OOE+OO 

37 

1. 

2.00E-06 

O.OOE+00 

O.OOE+OO 

38 

1. 

3.00E-06 

O.OOE+00 

O.OOE+OO 

39 

8. 

3.30E-06 

-4.00E+00 

-5.03E+03 

40 

7. 

2.00E-31 

-4.40E+00 

O.OOE+OO 

41 

1. 

l.OOE-09 

O.OOE+00 

O.OOE+OO 

42 

8. 

l.OOE-08 

-5.40E+00 

-2.10E+03 

43 

7. 

2. JOE-31 

-4.40E+00 

O.OOE+OO 

44 

1. 

l,00E-09 

O.OOE+00 

O.OOE+OO 

45 

1. 

l.OOE-09 

O.OOE+00 

O.OOE+OO 

46 

7. 

2.00E-31 

-4.40E+00 

O.OOE+OO 

47 

1. 

l.OOE-09 

O.OOE+OO 

O.OOE+OO 

48 

1. 

l.OOE-09 

O.OOE+00 

O.OOE+OO 

49 

8. 

1.20E-01 

-4.00E+00 

-8.80E+03 

50 

8. 

2.00E-05 

-3.30E+00 

-5.50E+03 

51 

8. 

5.00E-08 

-5.40E+00 

-2.10E+03 

52 

8. 

2.00E-04 

-3.30E+00 

-5.50E+03 

53 

8. 

5.00E-07 

-5.40E+00 

-2.10E+03 

54 

8. 

2.00E-03 

-3.30E+00 

-5.50E+03 

55 

8. 

l.OOE-02 

-5.70E+00 

-6.80E+03 

56 

8. 

l.OOE-02 

-5.70E+00 

-8.10E+03 

57 

1. 

5.00E-10 

O.OOE+OO 

O.OOE+OO 

58 

8. 

1.40E-29 

-l.OOE+00 

-6.00E+02 

59 

1. 

l.OOE-31 

O.OOE+OO 

O.OOE+OO 

60 

1. 

1.50E-10 

O.OOE+OO 

O.OOE+OO 

REACTIONS 

H''’.H200H  +  E  =  H2O  +  H2O,  1 

H'^.H20  +  E  =  H2  +  oh,  1 

H'''.(H20)2  +  E  =  H2O  +  H2  +  OH,  1 

H'^.(H20)3  +  E  =  oh  +  XX,  1 

H'^.(H20)4  +  E  =  OH  +  XX,  1 

H'^.H20.N2  +  E  =  H  +  XX,  1 

H'^.H20.C02  +  E  =  H  +  XX,  1 

H'^.(H20)2C02  +  E  =  H  +  XX,  1 

N0''‘.C02  +  E  =  no  +  CO2,  1 

N0'^.H20  +  I  =  NO  +  H2O,  1 

N0'*’.(H20)2  +  E  =  no  +  XX,  1 

N0'*’.(H20)3  +  E  =  no  +  XX,  1 

Oit"*"  +  M  =  02^  +  O2  +  M,  8 

NO'*'  +  N2  +  N2  =  N0'''.N2  +  N2,  7 

N0‘*’.N2  +  CO2  =  N0‘^.C02  +  N2,  1 

N0'''.N2  +  N2  =  no''"  +  N2  +  N2,  8 

N0'''.H20  +  N2  +  N2  =  N0'*’.H20.N2  +  N2,  7 

N0'^.H20,N2  +  CO2  =  N0'''.H20C02  +  N2,  1 

N0'^.H20.C02  +  H2O  =  NO‘'’.(H20)2  +  CO2,  1 

NO'*'.(H20)2  +  N2  +  N2=  N0'^.(H20)2.N2  +  N2 ,  7 

N0'''.(H20)2.N2  +  C02=  NO"'’.  (H20)2.C02  +  N2,  1 

N0'*'.(H20)2.C02  +  H2O  =  N0'''.(H20)3  +  CO2,  1 

H'’'.(H20)4  +  M  =  H'''.(H20)3  +  H2O  +  M,  8 

N0‘^.C02  +  M  =  NO"''  +  CO2  +  M,  8 

N0'^.H20.N2  +  M  =  N0'^.H20  +  N2  +  M,  8 

N0'''.H20.C02  M  =  N0'^.H20  +  CO2  +  M,  8 

N0'*'.(H2  0)2.N2  +  M  =  N0'^.(H20)2  +  N2  +  M,  8 

N0'''.(H20)2.C02+  M  =  N0'^.(H20)2  +  CO2  +  M,  8 

N0'^.(H20)3  +  M  =  N0‘*'.(H20)2  +  H2O  +  M,  8 

N0‘^.(H20)2  +  M  =  N0'^.H20  +  H2O  +  M,  8 

O2  +03=03  +  O2  ,  1 

O2  +  E  +  O2  =  O2  +  O2  ,  8 

O2  +  E  +  N2  =  O2  +  N2 ,  1 

O2  +  0  =  O3  +  E ,  1 
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TYPE 

PARI 

PAR2 

PAR3 

REACTION 

61 

1. 

1.50E-10 

O.OOE+OO 

0.001+00 

02~  +  0  =  0“  +  02,  1 

62 

1. 

2.00E-10 

O.OOE+00 

0. 001. +00 

O2  +  O2  ( ^  A )  =  O2  +  O2  + 

E, 

63 

7. 

2.50E-30 

-2.50E+00 

o.ooi;+oo 

O2  +  02  +  02  =  01+  +  O2 

,  7 

64 

1. 

4.30E-10 

O.OOE+OO 

O.OOE+OO 

04“  +  CO2  =  CO4"  +  02,1 

65 

1. 

1.30E-10 

O.OOE+OO 

O.OO.E+OO 

C04“  +03=  03~  +  CO2  +  1 

O2, 

66 

1. 

4.80E-11 

O.OOE+OO 

O.OOi+OO 

004“  +  NO  =  N03“  +  CO2, 

1 

67 

1. 

1.50E-10 

O.OOE+OO 

0.005+00 

CO4"  +  0  =  C03~  +  O2,  1 

68 

1. 

3.20E-10 

O.OOE+OO 

0.00 '5+00 

O3  +  0  =  O2  +  O2 ,  1 

69 

1. 

5.50E-10 

O.OOE+OO 

O.OOE+OO 

O3  +  CO2  =  CO3  +  O2 ,  1 

70 

1. 

l.lOE-10 

O.OOE+OO 

O.OOE+OO 

003“  +  0  =  02“  +  CO2,  1 

71 

1. 

l.lOE-11 

O.OOE+OO 

O.OOE+OO 

C03"  +  NO  =  N02''  +  CO2, 

1 

72 

1. 

2.00E-10 

O.OOE+OO 

O.OOE+OO 

CO3  +  NO2  =  NO3  +  CO2 , 

1 

73 

1. 

9.00E-11 

O.OOE+OO 

O.OOE+OO 

NO2  +03=  NO3  +  O2 ,  1 

74 

7. 

5.00E-30 

-1.90E+00 

O.OOE+OO 

0  +02+02=03  +  O2 , 

7 

75 

1. 

2.00E-10 

O.OOE+OO 

O.OOE+OO 

0~  +  0  =  O2  +  E,  1 

76 

1. 

3.00E-10 

O.OOE+OO 

O.OOE+OO 

0  +02(^4)  =03  +E,  1 

77 

1. 

4.40E-10 

O.OOE+OO 

O.OOE+OO 

0"  +  O3  =  O3"  +  0,  1 

78 

1. 

4.40E-10 

O.OOE+OO 

O.OOE+OO 

0  +03=02  +  O2 ,  1 

79 

1. 

6.00E-08 

O.OOE+OO 

O.OOE+OO 

02'*'  +  [NSUM]  =  XX,  1 

80 

1. 

6.00E-08 

O.OOE+OO 

O.OOE+OO 

04'*'  +  [NSUM]  =  XX,  1 

81 

1. 

6.00E-08 

O.OOE+OO 

O.OOE+OO 

O2  +  H2O  +  [NSUM]  =  XX, 

1 

82 

1. 

6.00E-08 

O.OOE+OO 

O.OOE+OO 

H'^.H20.0H  +  [NSUM]  =  XX, 

1 

83 

1. 

6.00E-08 

O.OOE+OO 

O.OOE+OO 

H‘*'.H20  +  [NSUM]  =  XX,  1 

84 

1. 

6.00E-08 

O.OOE+OO 

O.OOE+OO 

H'''.H20.N2  +  [NSUM]  =  XX, 

1 

85 

1. 

6.00E-08 

O.OOE+OO 

O.OOE+OO 

H‘^.H20.C02  +  [NSUM]  =  XX 

.  1 

86 

1. 

6.00E-08 

O.OOE+OO 

O.OOE+OO 

H.(H20)2'''  +  [NSUM]  =  XX, 

1 

87 

1. 

6.00E-08 

O.OOE+OO 

O.OOE+OO 

H.(H20)2.C02  +  [NSUM]  = 

XX, 

88 

1. 

6.00E-08 

O.OOE+OO 

O.OOE+OO 

(820)3  +  [NSUM]  =  XX, 

1 

89 

1. 

6.00E-08 

O.OOE+OO 

O.OOE+OO 

H'*’. (820)3  .CO2  +  [NSUM]  = 

XX 

90 

1. 

6.00E-08 

O.OOE+OO 

O.OOE+OO 

8'*'. (820)4  +  [NSUM]  =  XX, 

1 

91 

1. 

6.00E-08 

O.OOE+OO 

O.OOE+OO 

NO"^  +  [NSUM]  =  XX,  1 

92 

1. 

6.00E-08 

O.OOE+OO 

O.OOE+OO 

N0'^.N2  +  [NSUM]  =  XX,  1 

93 

1. 

6.00E-08 

O.OOE+OO 

O.OOE+OO 

N0'''.C02  +  [NSUM]  =  XX,  1 

94 

1 . 

6.00E-08 

O.OOE+OO 

O.OOE+OO 

02~  +  [PSUM]  =  XX,  1 
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TYPE 

PARI 

PAR2 

PAR3 

REACTION 

95 

1. 

6.00E-08 

O.OOE+00 

O.OOE+00 

0"  +  [PSUM]  =  XX,  1 

96 

1. 

6.00E-08 

O.OOE+00 

O.OOE+00 

04~  +  [PSUM]  =  XX,  1 

97 

1. 

6.00E-08 

O.OOE+00 

O.OOE+OO 

03 ~  +  [PSUM]  =  XX,  1 

98 

1. 

6.00E-08 

O.OOE+00 

O.OOE+OO 

C03~  +  [PSUM]  =  XX,  1 

99 

1. 

6.00E-08 

O.OOE+00 

O.OOE+OO 

CO4"  +  [PSUM]  =  XX,  1 

100 

1. 

6.00E-08 

O.OOE+00 

O.OOE+OO 

N02”  +  [PSUM]  =  XX,  1 

101 

1. 

6.00E-08 

O.OOE+00 

O.OOE+OO 

NO3"  +  [PSUM]  =  XX,  1 

102 

1. 

6.00E-08 

O.OOE+00 

O.OOE+OO 

N0'^.H20  +  [NSUM]  =  XX,  1 

103 

1. 

6.00E-08 

O.OOE+00 

O.OOE+OO 

N0'^.(H20)2  +  [NSUM]  =  XX,  1 

lOA 

1. 

6.00E-08 

O.OOE+00 

O.OOE+OO 

N0'^.(H20)3  +  [NSUM]  =  XX,  1 

105 

1. 

6.00E-08 

O.OOE+00 

O.OOE+OO 

N0'^.H20.N2  +  [NSUM]  =  XX,  1 

106 

1. 

6.00E-08 

O.OOE+00 

O.OOE+OO 

N0''’.H20.C02  +  [NSUM]  =  XX,  1 

107 

1. 

6.00E-08 

O.OOE+00 

O.OOE+OO 

N0'^.(H20)2.N2  +  [NSUM]  =  XX, 

108 

1. 

6.00E-08 

O.OOE+00 

O.OOE+OO 

N0'^.(H20)2.C02  +  [NSUM]  =  XX 

Table  2 

List  of  Reactloas  for  the  Sodium  Species 


TYPE^ 

.  ARl 

PAR2 

PAR3 

REACTIONS 

1 . 

3.00E-10 

O.OOE+00 

O.OOE+OO 

Na  +  O3  =  NaO  +  O2 ,  1 

1. 

4.00E-11 

O.OOE+00 

O.OOE+OO 

NaO  +  0  =  Na  +  O2 ,  1 

2. 

8.00E-31 

-2.90E+02 

O.OOE+OO 

Na^02  +  M  =  Na02  +  M,  2 

1. 

l.OOE-13 

O.OOE+00 

O.OOE+OO 

Na02  +  0  =  NaO  +  O2 ,  1 

1. 

O.OOE+00 

O.OOE+00 

O.OOE+OO 

NaO  +  H2O  =  NaOH  +  OH,  1 

1. 

5.00E-11 

O.OOE+00 

O.OOE+OO 

NaO  +  03=  Na02  +  O2 ,  1 

1. 

l.GOE-14 

O.OOE+00 

O.OOE+OO 

NaO  +  H  =  Na  +  OH,  1 

1. 

l.OOE-11 

O.OOE+00 

O.OOE+OO 

Na02  +  OH  =  NaOH  +02,1 

1. 

l.OOE-13 

0,00E+00 

O.OOE+OO 

Na02  +  H  =  NaOH  +  0,  1 

1. 

1.40E-12 

O.OOE+00 

O.OOE+OO 

NaOH  +  H  =  Na  +  H2O,  1 

1. 

l.OOE-10 

O.OOE+00 

O.OOE+OO 

NaOH  +  oId  =  NaO  +  OH,  1 

1. 

1.40E-U 

O.OOE+OO 

O.OOE+OO 

Na  +  HO2  =  NaOH  +0,  1 

1. 

l.OOE-10 

O.OOE+00 

O.OOE+OO 

NaO  +  O3  =  Na  +  XX,  1 

L. 

l.OOE-11 

O.OOE+OO 

O.OOE+OO 

NaO  +  HO2  =  NaOH  +  XX,  1 

1. 

6.00E-10 

O.OOE+OO 

O.OOE+OO 

Na  +  02"^  =  Na"^  +  O2,  1 

1. 

8.00E-11 

O.OOE+OO 

O.OOE+OO 

Na  +  NO"*"  =  Na"^  +  NO,  1 

7. 

2.00E-31 

-1.50E+00 

O.OOE+OO 

Na"*”  +  N2  +  M  =  Na^.N2  + 

M, 

7 

8. 

9.00E-07 

-2.50E+00 

-4.75E+03 

Na‘^.N2  +  M  =  Na'^  +  N2  + 

N2 

,  8 

7. 

2.00E-29 

-1.50E+00 

O.COE+00 

Na"*"  +  CO2  +  M  =  Na‘'’.C02 

+ 

M,  7 

8. 

6 .OOE-06 

-2.50E+00 

-6.25E+03 

Na'''.C02  +  M  =  Na'*'  +  CO2 

+ 

M,  8 

7. 

5.00E-32 

-1.50E+00 

O.OOE+OO 

Na^  +  O2  +  M  =  Na02'*^  +  M 

1, 

7 

8. 

8.00E-09 

-2.50E+00 

-2.85E+03 

Na02'*'  +  M  =  Na"*”  +  O2  +  N 

1, 

8 

7. 

l.OOE-28 

-1.50E+U0 

O.OOE+OO 

Na^  +  H2O  +  M  =  Na^.H20 

+ 

M,  7 

1 . 

1 .OOE-09 

O.OOE+OO 

O.OOE+OO 

Na^.N2  +  CO2  =  Na^.C02  + 

•  N2,  1 

1 . 

l.OOE-09 

O.OOE+OO 

O.OOE+OO 

Na'^.C02  +  H2O  =  Na'''.H20 

+ 

CO2, 

7. 

1  .OOE-06 

-5.00E-01 

0.(!0E+00 

Na'''.C02  +  E  =  Na  +  CO2, 

7 

The  parameterization  and  values  for  the  reaction  constants  are  taken  from 
Swider  and  Narcisi  (1984).  The  reaction  rates  are  calculated  as  follows: 
Type  1:  k  =  PARI 

Type  2:  k  =  PARl*exp(-PAR2yT) 

Type  7:  k  =  PARl*( 

Type  8:  k  =  PARl*(  30U/T)~^''‘^^*exp( PAR3/T) . 


TYPE 

PARI 

PAR2 

PAR3 

REACTIONS 

27 

7. 

l.OOE-06 

-5.00E-01 

O.OOE+00 

Na'''.H20  +  E  =  Na  +  H2O 

28 

7. 

l.OOE-06 

-5.00E-01 

O.OOE+00 

Na'*'.N2  +  E  =  Na  +  N2, 

29 

7. 

l.OOE-06 

-5.00E-01 

O.OOE+00 

Na02''’  +  E  =  Na02  +  ,  7 

30 

1. 

l.OOE-10 

O.OOE+00 

O.OOE+00 

Na02'*’  +  0  =  NaO'*'  +  O2 , 

31 

1. 

5.00E-10 

O.OOE+00 

O.OOE+00 

NaO'*'  +  0  =  Na"*"  +  O2,  1 

32 

1. 

l.OOE-09 

O.OOE+00 

O.OOE+00 

NaO  +  02+  =  NaO"*"  +  O2 , 

33 

1. 

l.OOE-09 

O.OOE+00 

O.OOE+00 

NaO  +  NO"*"  =  NaO"'’  +  NO, 

34 

1. 

4.00E-07 

O.OOE+00 

O.OOE+00 

NaO'*’  +  E  =  Na  +  0,  1 

35 

1. 

6.00E-08 

O.OOE+00 

O.OOE+00 

Na"*"  +  [NSUM]  =  Na  +  XX 

36 

1. 

6.00E-08 

O.OOE+00 

O.OOE+00 

Na  +  N2  +  [NSUM]  =  Na  ■ 

37 

1. 

6.00E-08 

O.OOE+00 

O.OOE+00 

Na  +  CO2  +  [NSUM]  =  Na 

38 

1. 

6.00E-08 

O.OOE+00 

O.OOE+OO 

Na  +  H2O  +  [NSUM]  =  Na 

39 

1. 

6.00E-08 

O.OOE+00 

O.OOE+OO 

NaO'*'  +  [NSUM]  =  NaO  +  : 

40 

1. 

6.00E-08 

O.OOE+00 

O.OOE+OO 

Na02''‘  +  [NSUM]  =  Na  + 
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Table  3 

Calculated  Local  Topside  Scale  Heights  (km) 


1 

Altitude  1 

Scale  Height! 

(km) 

(km)  j 

No  ions 

K*  =  7  X  10"31tt  K*  = 

2  X  10*31 

K*  =  4  X  10"32 

1 

92  1 

9.2 

8.13 

7.84 

6.85 

94  1 

6.72 

6.52 

6.11 

5.23 

96  1 

6.13 

5.74 

5.12 

4.31 

98  1 

5.85 

4.99 

4.27 

3.31 

100  1 
_ L 

5.73 

4.26 

3.60 

2.48 

tLocal  scale  height  is  obtained  from  the  gradient  of  Na  profile. 

ttNa^  +  N2  +  M  >  Na^  ->■  N2  +  M  with  reaction  rate  Kg  =  K*  [■^— )  .  Units  for 

K*  are  cm®  s  ^  . 


Experimental  Results 

Gibson  and  Sandford  (1970):  4  km  (98  km),  3  km  (100  km),  2  km  (102  km). 

Megie  and  Blamont  (1977):  2-5  km  average  scale  height  above  peak. 

Siraonich  et  al.  (1979):  2.5-5  km  at  100  km  (scale  height  changes  during 
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Figure  Captions 


Figure  1.  Schematic  diagram  for  the  chemical  scheme  of  the  positive  ions  in 
the  D  and  lower  E  regions*  The  icns  are  grouped  into  two  families 
with  02^  and  NO^  being  the  primary  ions  in  each  group. 

Figure  2.  Schematic  diagram  for  the  chemical  scheme  of  the  negative  ions  in 
the  D  and  lower  E  regions. 

Figure  3.  The  production  rates  for  positive  ions  for  daytime  and  nighttime 

conditions:  A)  Day  (45®  Solar  Zenith  Angle  (SZA)),  altitude  range 
90-120  km);  B)  Day  (45®  SZA),  altitude  range  60-90  km;  c)  Night 
(180®  SZA),  altitude  range  60-120  km.  An  altitude-independent 
production  rate  of  10~3  cm^  s~l  due  to  gallactic  x-rays  is  also 
included  in  the  model.  The  curves  are:  a)  x-ray,  b)  Lyman-B,  c) 
Lyman-a ,  d)  energetic  electrons,  e)  ionization  of  02(^A),  and  f) 
cosmic  ra>s. 

Figure  4.  Flow  chart  for  the  diurnal  model  of  neutral  and  ionic  sodium 
species. 

Figure  5.  Calculated  diurnal  variation  of  electron  density.  Profiles  are 
calculated  at  (a)  86  km,  (b)  90  km,  (c)  94  km,  and  (d)  100  km. 

Figure  6.  Calculated  diurnal  vari  ition  of  the  densities  of  NO^  and  02^. 

Profiles  are  calculated  at  (a)  86  km,  (b)  90  km,  (c)  94  km,  and 
(d)  100  km. 

Figure  7.  Calculated  total  electron  and  negative  ion  densities  for  daytime 

and  nighttime.  The  positive  ion  densities  are  equal  to  the  sum  of 
the  electron  and  the  negative  ion  densities.  Curves  (a)  and  (b) 
are  the  calculated  alti  ;ude  profiles  for  electrons  for  day  and 
night,  respectively.  Tie  corresponding  day  and  night  profiles  for 
the  negative  ions  are  given  by  curves  (c)  and  (d).  Curves  (e)  and 
(f)  are  the  profiles  for  electron  and  negative  ions  calculated  for 
nighttime  including  the  production  of  ions  from  energetic 
electrons . 

Figure  8.  Calculated  altitude  profiles  of  negative  ions  for  day  (A)  and 

night  (B).  The  curves  are  02~  ( - ),  CDs”  ( - ),  C0n~  (***),  NOa" 

(-•-•-),  N02"  (-0-0-),  0“  (-A-4-). 

Figure  9.  Calculated  altitude  profiles  of  02^  ("a"  for  day,  "b"  for  night) 
and  NO’*'  ("c"  for  day  and  "d"  for  night). 

Figure  10.  Calculated  day  (A)  and  night  (B)  altitude  profiles  of  the  dominant 

positive  ions  below  100  km.  The  curves  are  H'*’.(H20)4  ( - ), 

N0'*^.(H20)  ( - ),  N0^.(H20)2  (•••)»  N0^.(C02)  (-•-•-),  H^(H20)3 

(-0-0-). 

Figure  11.  Schematic  diagram  for  the  chemical  scheme  of  sodium  neutral  and 
ionic  species. 
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Figure  12.  Calculated  diurnal  variation  of  the  Na'*’  density  at  (a)  90  km,  (b) 
94  km,  and  (c)  100  km. 

Figure  13.  Calculated  altitude  profiles  for  Na  and  Na'*’  at  noon. 

Figure  14.  Calculated  altitude  profiles  of  atomic  sodium  using  different 

rates  for  reaction  (4-4).  The  curves  correspond  to  values  of  K* 

2  X  10"31  cvfi  s"l  ( _ ),  K*  »  7  X  10“31  cm^  s"l  ( ),  K*  -  4  x 

10~32  cm^  s~2  (-X-),  and  the  case  where  no  Na"*"  is  present  (-o-o-) 

Figure  15.  Calculated  profiles  of  the  Na  layer,  assuming  vertical  ion  drift 
velocities  of  0  ( - ),  -4  (-A-A-),  -10  (-o-o-),  and  -20  (-x-x-) 
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Figure  3.  The  production  rates  for 
positive  ions  for  daytime  and 
nighttime  conditions:  A)  Day 
(45®  Solar  Zenith  Angle  (SZA)), 
altitude  range  90-120  km);  B) 
Day  (45°  SZA),  altitude  range 
60-90  km;  c)  Night  (180®  SZA), 
altitude  range  60-120  km.  An 
altitude-independent  production 
rate  of  10”^  cm^  s”^  due  to 
gallactic  x-rays  is  also  in¬ 
cluded  in  the  model.  The 
curves  are:  a)  x-ray,  b) 
Lyman-6 ,  c)  Lymah-o ,  d)  ener¬ 
getic  electrons,  e)  ionization 
of  02(^4),  and  f)  cosmic  rays. 
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Calculated  diurnal  variation  of  the  densirlci  of  NO'*’  and  Oz*- 
Profiles  are  calculated  at  (a)  86  km,  (b)  90  kra,  (c)  9A  kra,  and 
(d)  100  kni> 
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Calculated  total  electron  and  negative  ion  densities  for  daytime 
and  nighttime.  The  positive  ion  densities  are  equal  to  the  sum  of 
the  electron  and  the  negative  ion  densities.  Curves  (a)  and  (b) 
are  the  calculated  altitude  profiles  for  electrons  for  day  and 
night,  respectively.  The  corresponding  day  and  night  profiles  for 
the  negative  ions  are  given  by  curves  (c)  and  (d).  Curves  (e)  and 
(f)  are  the  profiles  for  electron  and  negative  ions  calculated  for 
nighttime  including  the  production  of  ions  from  energetic 
electrons . 
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Figure  9. 


Calculated  altitude  profiles  of  O2'*’  ("a”  for  day,  "b"  for  night) 
and  NO''’  ("c”  for  day  and  "d"  for  night). 
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Figure  10. 


Calculated  day  (A)  and  night  (B)  altitude  profiles  of  the  dominant 
positive  ions  below  100  km.  The  curves  are  H  ^ 

NO-'.CHaO)  ( - ),  N0‘".(H20)2  (•••).  N0^.(C02)  ( - ),  H  (H20)3 

(-0-0-). 
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Figure  12.  Calculated  diurnal  variation  of  the  Na"*"  density  at  (a)  90  km,  (b) 
94  km,  and  (c)  100  km. 
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Figure  14.  Calculated  altitude  profiles  of  atomic  sodium  using  different 

rates  for  reaction  (4-4).  The  curves  correspond  to  values  of  K* 

2  X  10~31  cn^  s~l  ( - ),  K*  -  7  X  lO'^l  cm^  s"^  ( - ),  K*  -  4  x 

10“^2  cm®  s”2  (_x-)^  and  the  case  where  no  Na"'"  Is  present  (-o-o-) 
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PART  B 


Impact  of  Neutral  Sodium  Chemistry  on  the  Stratospheric  Chlorine  Cycle 


Recent  theoretical  and  experimental  developments  in  sodium  chemistry  have 
created  much  interest  in  its  role  in  stratospheric  chlorine  chemistry.  Murad, 
Swider,  and  Benson  (1981)  first  suggested  that  the  reaction 


*R19: 


NaOH  +  HCl 


NaCl  +  H^O 


would  provide  a  favored  pathway  for  removal  of  NaOH  in  competition  with  NaOH 
photolysis  if  >  10”^^  cm^  s~^ .  If  NaCl  does  not  react  to  release  the  Cl 
atoms,  R19  could  be  an  Important  sink  for  stratospheric  chlorine  species  if 
the  atmospheric  concentration  of  total  sodium  (NaX)  is  sufficiently  high. 

On  the  other  hand,  Rowland  and  Rogers  (1982)  suggested  that  NaCl  could  be 
rapidly  photolysed  based  on  the  assumption  of  unity  quantum  yield  and  the 
absorption  cross  sections  measured  at  high  temperatures.  Reaction  R19, 
followed  by 


R20:  NaCl  +  hv 


Na  +  Cl 


could  then  provide  an  additional  mechanism  for  converting  HCl  to  Cl,  similar 
to  the  well-known  source 


OH  +  HCl  — ^  HjO  +  Cl 


at  altitudes  below  50  km. 

Subsequent  measurements  by  Silver  et  al.  (1984)  have  obtained  values  of 
2.5  X  10“^°  cra^  s”^  for  reaction  R19.  This  reaction  will  affect  the  parti¬ 
tioning  of  the  chlorine  species  in  the  upper  stratosphere  when 


♦Reaction  numbers  R  and  rates  k  ,  J  ,  refer  to  the  reactions  in  Table  1. 
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Jgo  [NaCl]  ~  k^[OH] [HCl] 


(2) 


If  Ri9  is  the  only  mechanism  for  forming  NaCl,  (2)  reduces  to 


kj  [NaOH]  ~  k^IOH] 


which  translates  into 


[NaOH]  ~  IC/*  cm“3 


for  typical  stratospheric  concentration  of  10^  cm~3  for  OH. 

The  density  of  sodium  species  NoX  (Na,  NaO,  Na02»  NaOH,  NaCl)  in  the 
lower  stratosphere  is  ~  10^  cm~3 ,  if  one  assumes  that  sodium  diffuses  down 
from  the  mesosphere  and  is  removed  in  the  troposphere  with  a  time  constant  of 
ten  days  (Liu  and  Reid,  1979;  Sze  et  al.,  1982).  In  the  stratosphere,  most  of 
the  NaX  is  expected  to  be  in  the  form  of  Na02  and  NaOH.  Past  attempts  in  as¬ 
sessing  the  impact  of  Na  on  stratospheric  chlorine  chemistry  were  hampered  by 
the  lack  of  kinetic  data  to  determine  the  NaOH  concentration  in  the  stratosphere. 

Recent  measurements  reported  by  Silver  and  Kolb  (1985)  indicated  a  rate 
of  2  X  10”^ 0  cm^  s~^  for  the  reaction 


R2 1 ;  Na02  +  HCl  ->■  NaCl  +  HO2  . 


Since  both  NaOH  and  Na02  react  with  HCl  at  about  the  same  rate,  the  partition¬ 
ing  between  the  two  reservoirs  is  no  longer  crucial.  Indeed,  the  concentra¬ 
tion  of  NaCl  in  the  stratosphere  can  be  estimated  from  the  total  sodium  con¬ 
centrations  via  the  following  set  of  algebraic  equations: 


J20[NaCl]  =  k*[HCl][NaR] 
[NaCl]  +  [NaR]  S  [NaX] 


(5a) 


(5b) 


where  [NaR]  is  the  sum  of  [NaOH]  and  tNa02],  and  k*  =  2  x  cm^  s”l  is  the 

reaction  rate  constant  for  Na02  and  NaOH  with  HCl.  Equation  (5a)  follows  from 
the  photochemical  equilibrium  assumption  for  NaCl  while  equation  (5b)  follows 
from  the  closure  relationships  assuming  that  [Na]/[NaX]  and  [NaOl/[NaX]  <  1. 
From  equations  (5),  it  follows  that 


[NaCl]  = 


k*[HCl] [NaX] 
Jon  +  k*[HCl' 


(6) 


In  the  stratosphere  (Z  <  50  km)  with  [HCl]  ~  10®  cm~® ,  k*(HCl]  ~  10”^  s'^ .  As 

long  as  J  <  k*[HCl],  the  sodium  species  will  be  in  the  form  of  NaCl.  Thus, 

the  effectiveness  of  R20  in  recycling  HCl  to  Cl  will  be  linear  in  J20  and 

[NaX]. 

Our  present  knowledge  of  the  reactions  controlling  sodium  chemistry  is 
summarized  in  T^ble  1.  The  photolysis  rates  for  NaCl  and  NaOH  (reactions  R20 
and  R4),  estimated  by  Rowland  and  Rogers  (1982)  and  Rowland  and  Makide  (1982) 
can  be  considered  only  upper  bounds  to  the  actual  rates.  Rates  for  other 
reactions  are  usually  estimated  from  analogy  to  similar  reactions  involving  H 
instead  of  Na.  Observations  of  the  properties  of  the  mesospheric  sodium  layer 
and  the  sodium  D-llne  nlghtglow  are  used  to  constrain  the  rate  for  reactions 
R2,  R5,  and  R6  (Bates  and  Ojha,  1980;  Rodriguez,  Ko,  and  Sze,  1984). 

We  will  now  present  some  quantitative  results  using  our  stratospheric 
model.  The  value  adopted  for  J20  in  the  model  is  shown  in  Figure  1  and  is 
similar  to  those  suggested  by  Rowland  and  Rogers  (1980).  Figure  2  shows  the 
calculated  profiles  for  Na02 .  NaOH,  and  NaCl.  Figure  3  compares  the  rate  of 
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recycling  HCl  to  Cl  by  photolysis  of  NaCl  with  the  reaction  of  OH  +  HCl.  The 
calculated  conversion  rate  is  linear  in  the  J20  and  [NaX]  in  the  lower 
stratosphere,  and  the  mechanism  could  dominate  the  OH  +  HCl  reaction  if  J20  is 
sufficiently  fast  (i.e.,  J20  ~  10”^  s  ^). 

Note  that  the  above  calculation  should  be  considered  as  the  upper  limits 
since  heterogeneous  reactions  involving  NaCl,  or  dimerization  (Lamb  and  Ben¬ 
son,  1985),  may  prevent  it  from  photolysing.  In  addition,  reaction  of  NaOH 
with  other  species  may  prevent  formation  of  NaCl.  One  possibility  was  sug¬ 
gested  by  Murad  and  Swider  (1979)  where  NaOH  reacts  with  CO2  to  form  NaHCOs : 

NaOH  +  CO2  NaHC03  . 

Estimated  reaction  rate  is  ~  10~32  ^m^  ,  giving  a  loss  frequency  for  NaOH 

of  about  2  X  10~2  s'^ ,  .475  s"! ,  and  10  s~^ ,  at  40  km,  30  km,  and  20  km 
respectively.  If  the  reaction  indeed  proceeds  at  this  rate,  all  the  sodium 
species  will  be  in  the  form  of  NaHCOs . 

Although  there  are  large  uncertainties  in  the  sodium  chemistry,  the  fol¬ 
lowing  statements  can  be  established.  Given  the  stratospheric  chlorine  con¬ 
centration  of  about  10®  to  10®  cm~® ,  the  present  estimated  concentration  of 
10®  cm~®  for  NaCl  will  not  constitute  an  important  reservoir  for  chlorine. 

Any  impact  on  the  chlorine  cycle  must  then  involve  a  catalytic  mechanism.  Our 
calculation  shows  that  if  gaseous  NaCl  is  the  major  sodium  compound  and  J20  > 
10~®  s  photolysis  of  NaCl  will  be  an  Important  pathway  for  cycling  HCl  to 

Cl  in  the  stratosphere.  However,  the  calculation  did  not  take  into  considera¬ 
tion  the  possibility  of  the  formation  of  other  sodium  compounds  such  as 
NaHCOa ,  (NaCl)  ,  or  incorporation  of  NaCl  into  aerosols. 


Table  1:  Reaction  Scheme  and  Rate  Constants 


Rate^ 


Reaction 


1 . 

Na  +  O3 

-► 

NaO  +  O2 

3.5(-10) 

[h] 

2. 

NaO  +  0 

Na  (2p,  2s)  +  O2 

4(-ll) 

[c] 

3. 

Na  +  O2  +  N2 

Na02  +  N2 

1.9  X  10"30(T/300)“1*1 

[d] 

4. 

NaOH  +  hv 

Na  +  OH 

2(-3) 

[e] 

5. 

Na02  +  hv 

-► 

Na  +  O2 

l(-4) 

[c] 

6. 

Na02  +  0 

4- 

NaO  +  O2 

1(-13) 

[c] 

7. 

NaO  +  H2O 

> 

NaOH  +  OH 

2(-10) 

[i] 

8. 

NaO  +  O3 

> 

Na02  +  O2 

1.8(-10) 

[il 

9. 

NaO  +  H 

-► 

Na  +  OH 

1(-14) 

[g] 

10. 

Na02  +  OH 

-► 

NaOH  +  O2 

K-ll) 

[f] 

11. 

Na02  +  H 

> 

NaOH  +  0 

1(-13) 

[g] 

12. 

NaOH  +  H 

> 

Na  +  H2O 

1.4(-12) 

[f] 

13. 

NaOH  +  (Id) 

-► 

NaO  +  OH 

K-lO) 

[g] 

14. 

Na  +  HO2 

-f 

NaOH  +  0 

K-iO) 

[f] 

15. 

NaO  +  O3 

-f 

Na  +  202 

K-lO) 

[1] 

16. 

NaO  +  HO2 

*► 

NaOH  +  O2 

K-ll) 

[f] 

17. 

NaO  +  H2 

-► 

Na  +  H2O 

i(-ll) 

[i] 

18. 

NaO  +  H2 

-► 

NaOH  +  H 

2.6(-ll) 

[i] 

19. 

NaOH  +  HCl 

-► 

NaCl  +  H2O 

2.5(-10) 

[d] 

20. 

NaCl  +  hv 

> 

Na  +  Cl 

J  <  6(-2) 

[j] 

21. 

Na02  +  HCl 

-► 

NaCl  +  HO2 

2.5(-10) 

[K] 

‘Read 

3(-10)  as  3 

X 

10~1*^.  Units  are  cm^  s~^ , 

except  cm^  s~l  for  (3)  and 

1  s~ 

for  (4),  (5),  and  (15).  [b]  Kolb  and  Elgin  (1976);  [c]  Sze  et  al.  (1982); 

[d]  Silver  et  al.  (1984);  [e]  Rowland  and  Makide  (1982);  [f]  Liu  and  Reid 
(1979);  [g]  Kirchoff  and  Clemesha  (1983);  [h]  Zahniser  et  al.  (1985);  [i]  Age 
and  Howard  (1985);  [j]  Rowland  and  Roger  (1982);  [k]  Silver  and  Kolb  (1985). 
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Figure  Captions 


Figure  1.  Photolysis  rate  of  NaCl  used  in  the  model  calculation.  The  values 
are  taken  from  Rowland  and  Rogers  (1982). 

Figure  2.  Calculated  concentrations  of  the  sotlium  species  in  the  strato¬ 
sphere.  The  concentrations  of  Na  and  NaO  are  less  than  10  cra”3  in 
the  lower  stratosphere. 

Figure  3.  Comparison  of  the  conversion  rate  of  HCl  to  Cl  by  NaCl  photolysis 
( - )  and  reaction  OH  +  HCl  ( - ). 
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